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Abstract 

The kink signature of charged tracks is predicted in some SUSY models, and it is 
very characteristic signal at collider experiments. We study the kink signature at LHC 
using two models, SUSY models with a gravitino LSP and a stau NLSP, and R-parity 
violating SUSY models with a stau (N)LSP. We find that a large number of kink events 
can be discovered in a wide range of the SUSY parameters, when the decay length is 
0(10 - 10 5 ) mm. Model discrimination by identifying the daughter particles of the kink 
tracks is also discussed. 



1 Introduction 



In most of supersymmetric (SUSY) standard models, it is assumed that the R-parity is exactly 
conserved, and the lightest SUSY particle (LSP) is the lightest neutralino. The conservation 
of the R-parity makes the LSP neutralino stable and promotes it to a viable candidate for 
the dark matter. A characteristic signature for such a scenario at the Large Hadron Collider 
(LHC) is a large missing transverse momentum. 

However, some SUSY models provide not such a signature but more exotic ones. A 
popular example is the SUSY models with a gravitino LSP, such as in gauge-mediated SUSY 
breaking models O. In this case, the next-to-lightest SUSY particle (NLSP) becomes long- 
lived due to the extremely weak interaction of the gravitino. Another example is the R-parity 
violation 0. Recently some of the authors studied cosmological constraints on the R-parity 
violation taking into account the effect of lepton flavor violations, and showed that the baryon 
asymmetry of the universe favors very tiny R-parity violations ||. This means that the LSP 
(or NLSP, if the LSP is the gravitino) decays into Standard Model particles with a long decay 
length. 

The long-lived superparticles in such scenarios provide various kinds of interesting 
signatures at the LHC The signatures crucially depend on the decay length cr and the 

electric charge of the long-lived superparticle X: 

(i) Stable signature: For cr > L, where L ~ 0(m) is a typical detector size, most of the 
produced X particles pass through the detector before decaying. If X is electrically 
neutral, the signature is a missing transverse momentum, which is similar to the 
conventional SUSY scenario. However, when X is electrically charged, as in the case of 
long-lived stau scenario, it provides an anomalous track of large ionization and/or low 
velocity Q. Furthermore, the charged track can be used to determine the masses of 
X [Q-[iIJ and other superparticles J|, 10, 11 1, as well as the spin and other properties of 



superparticles [12|. The lifetime of X may also be measured by looking for its delayed 
decay after being stopped |jl3|| . Recently CMS has reported the lower limit of the 
stable/stopped gluino mass, 398/370 GeV respectively, using the highly ionized track 



and empty bunch jl4, 15 1. ATLAS also reported lower limit on the stable gluino mass, 



562 GeV |§. 



(ii) In-flight decay: When the decay length satisfies 0(cm) < cr < NL, where N is the 
total number of X produced through the experiment, its in-flight decays can be seen 
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within the detectors g. [Tt] fHJ .Q In the case of neutralino NLSP decaying into gravitino 
LSP, non-pointing photons may be observed [17]. If X is charged, its in-flight decay can 



cause a disappearing track |L8|,|2C|] or a kinked charged track. The lifetime measurement 
may also be possible by studying the distribution of, or simply by counting the number 
of, these in-flight decay events [||, |l7|-|l9| . 



Note that these two cases are not always mutually exclusive; for L < cr < NL, both types 
of signatures can be studied simultaneously. 

In this paper, we will study the latter case (ii), in particular, the kinks of charged 
tracks which are produced when long-lived charged particles decay into another charged 
particle inside the detector. Specifically, we investigate the feasibility of probing the SUSY 
models by analyzing such kink tracks in the transition radiation tracker (TRT) of the ATLAS 
detector, which provides us continuous detection of charged tracks. As underlying models, 
we concentrate on the following two models; (a) SUSY models with a gravitino LSP and a 
stau NLSP, and (b) R-parity violating SUSY models with a stau (N)LSP. We show that, 
in a realistic setup, a large number of kink events can be discovered in a wide range of the 
SUSY parameters. We also discuss the possibility of model discrimination by identifying the 
daughter particles of the kink tracks. 

2 Kink Tracks 

A kink track signature and possible backgrounds are discussed here. When a particle decays 
in tracking detectors and produces one charged daughter particle, the track would bend 
abruptly, which forms a kink track. In the ATLAS detector, charged tracks are observed in 
the inner detector, which consists of the pixel detector, the semi-conductor tracker (SCT) 
and the TRT, located from the inside to outside of the beam axis respectively. In particular, 
the TRT, which is a continuous tracking detector, is suited to observe the kinks directly. 

The TRT detector [21] consists of the barrel part and the end-cap part. The barrel 
part covers the collision point cylindrically: \z\ < 712mm and 563mm < R < 1066mm, 
while the end-cap part caps both sides.^ The barrel part consists of three modules: the first 
module (the innermost one) sits at the range of 563 mm < R < 694 mm, the second one at 



1 The probability of the X particle decaying within the distance L is given by P(L) = 1 — exp(— L/P'jct), 
where /3 is the velocity of X and j = (1 — /3 2 ) _1//2 , and typically /3j = O(l). For cr 2> L, it is evaluated as 
P(L) ~ L/P^ct ~ L/ct, and hence the number of in-flight decay events is given by ~ (L/ct)N. 

2 Our coordinate system follows the ATLAS standard: z is along the beam axis whose origin is the nominal 
collision point, R is the distance from the beam axis, and and <j> is the azimuthal and the polar angles. 
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697 mm < R < 860 mm, and the third one (the outermost one) at 863 mm < R < 1066 mm. 
Each module contains straw tubes, which are aligned parallel to the beam axis and provide 
R- <j) information of the track position. The barrel TRT is composed of 73 layers of the straws 
in total. When a charged particle passes through the TRT, it hits a large number (typically 
~ 40) of the straws. In the analysis we use the barrel part of the TRT to measure the kink 
tracks. 

The kink track is identified by reconstructing both the tracks of the charged mother 
and daughter particles. Schematic pictures of the kink track are shown in Fig. [l|. The 
mother particle flies through the pixel and the SCT, and reaches the TRT. Consider that 
the mother particle decays in the TRT. The track terminates halfway, but is reconstructed 
very well by the information from the pixel detector and the SCT f20|| . A good momentum 
resolution |]20| , p2| is also expected, and it is improved by combining the information from the 
TRT Moreover, it may be possible to measure the lifetime and the mass of the mother 
particle (or their combination) by the information on the momentum, the time-of-flight and 



the decay length, depending on the models [TLS , 20 



The track of the mother particle terminates when the particle decays. As we will require 
the mother particle to be \r]\ < 0.63, where r] is the pseudorapidity, it would pass through 
the full volume of the 3rd module of the TRT if it does not decay. So we can recognize the 
termination of the track by identifying the absence of the TRT hits (less than noise level) 
after the decay point. Moreover, a number of hits will be observed in the 3rd module for the 
daughter particle(s). However, when the daughter particle has the same azimuthal angle as 
the mother, these two tracks can not be distinguished in the TRT. In order to reconstruct the 
tracks separately to observe a kink, the kink is required to be Acj) > 0.1, where the kink angle 
A(f> is defined as the difference of the azimuthal angles between the mother and daughter 
tracks. The charged daughter particle is also required not to escape from the end-cap of the 
TRT, which guarantees that the daughter particle passes through full volume of the TRT 
3rd module (about 15 hits are expected), and the daughter track can be reconstructed with 
reasonable efficiency. Using the simulations based on G4, in which all detector materials are 
taken into account, the reconstruction efficiency is estimated to be 70% for Pt > 20 GeV, 
and 60% even for ~ 10 GeV [|24|j. This reconstruction efficiency is taken into account in 
this study. 

We also need to comment that it is difficult or impossible to reconstruct the tracks of 



soft [20, 2^ or neutral daughter particles. Especially we can hardly reconstruct the kink for 
the case where the charged daughter particle is soft, which usually occurs when a neutral, 
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Figure 1: Schematic pictures of the decay which generates a kink track in the TRT. In the 
left panel, the tracks are projected onto the plane which includes the beam axis. In the right 
panel, they are projected onto the plane which is perpendicular to the beam axis. 



i.e. missing, daughter particle is almost degenerate with the mother particle. Even in such 
cases, the mother track can be identified as a disappearing track and we can measure the 



lifetime of the mother particle [20] . In this paper we consider models in which the decay 
products are much lighter than the mother particle. We can identify not only the mother but 
also the daughter particles. 

It is possible to identify the particle species of the daughter track from the information of 
various detectors. When the daughter particle is an electron, we obtain signals in the TRT 
and the electromagnetic (EM) calorimeter. When a muon, a characteristic signature arises in 
the muon spectrometer. When the mother particle decays into jets, they are identified by the 
EM and hadronic calorimeters. In these three cases, the momentum of the daughter can be 
determined. Lastly, when the daughter is a tau lepton, it decays right after the production 
and emits an electron, a muon, or hadrons. The ratio of the kink events with the e/^/jets 
final state would enable us to identify the daughter tau production statistically, which will be 



discussed in Sec. 3.3. Once the decay branching fraction of the mother particle is known, it 
helps us to discriminate the underlying models, such as the models with the R-parity violation 
or with a gravitino LSP. 

The masses of the mother particles and the invisible daughter can also be determined. 
Velocity of the mother particle can be measured with an accuracy of 0(10)% from pixel 
energy loss information |R|, and accuracy of 0(1)% from the timing information in the 
hadronic calorimeter and the muon detectors H^]. If the typical CT7 is longer than 0(1) m, 
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the mother particle's mass can be measured precisely. The transverse momentum Pr of 
the visible daughter particles can be measured in TRT and the kinematic edge of the Pp 
distribution carries the information of the invisible daughter particle's mass, with a typical 
resolution of 0(1) GeV. 

The possible background sources for the kink signature are as follows: 

1. In-flight decays of hadrons, such as tt^ and K^, have similar topology. However, it 
is considered that the hadronic background events are suppressed by requiring a high 
transverse momentum of the mother particle and a large kink angle. This is because 
a boosted particle is likely to emit the decay products forward. Consider a charged 
particle X decaying into a charged particle A and another (neutral) particle. The kink 
angle A(p (the difference of the azimuthal angles between X and A) has a kinematical 
bound 

• A^l 1 mxlp^l 1 m\/2 

I sm < s - m0x mA | p Lab| < sin6x mA P T (X) ' ( j 

where "Lab" denotes the laboratory frame, and "CM" does the rest frame of X. Here, 
mx and txia are the masses of X and A, respectively. The 3-momentum of the particle 
i is represented by pf in the frame P, while P r r(X) shows the transverse momentum 
of X in the laboratory frame. The angle Ox is the polar angle of X, that is, the angle 
between X and the beam axis. In the analysis, we consider the decay of the mother 
particle to be restricted in \rj\ < 0.63, which corresponds to 0.83 < sinOx < 1. In order 
to suppress the in-flight decay background, we impose Pt(X) > 100 GeV and A(p > 0.1. 
Using the simulation based on PGS, the contribution of in-flight decay is estimated to be 
negligibly small. 

2. Stable charged hadrons (n^, K^, p, etc.) interact hadronically with materials of the 
inner detectors and the direction of the particle is changed significantly. This process 
also makes kink signatures. However, background events with stable charged particles 
are also suppressed by the above requirement that the mother particle must have Pr > 
100 GeV. 

3. The unphysical tracks produced by noise hits in the pixel and SCT detectors can have 
Pr > 100 GeV, but these contributions can be suppressed dramatically by requiring the 
daughter track. 

Background rate is expected to be very low after all selections, which we will discuss later, 
are applied. 
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3 Models 



In this section we briefly introduce super symmetric models which can produce kink tracks at 
the LHC. We also discuss the model discrimination in Sec. 



3.1 Models with a gravitino LSP and a stau NLSP 

In super symmetric models with a gravitino LSP and a stau NLSP, under the R-parity 

conservation, the stau decays into a tau lepton and a gravitino and becomes long-lived. 
Its lifetime is given by 

/ 5 \ - 1 

/ m~ \ 

ct ~ 



2 

3/2 / 

2 



* 550mm ( ^ ) (H^/L) , (2) 
V200GeV/ VlkeV/ V ; 

where M p is the reduced Planck scale, is the stau mass, and 77i 3/ / 2 is the gravitino mass. 

As we shall see, the stau kink may be observed at the TRT for m 3/ / 2 ~ 0(0.01 - 1) keV. 

Interestingly, this corresponds to the most challenging gravitino mass range from cosmological 

point of view |q, p5|, which will be discussed in Sec. |j| The produced tau lepton decays 

promptly, either leptonically or hadronically. Then, the signature will be a kink track with 

the daughter particle being an electron, a muon, or a tau-jet. 

The stau can be long-lived also in the case of an axino LSP and a stau NLSP. In the 



hadronic, or KSVZ, axion models [26], the stau lifetime becomes much longer than the 



detector size [27|, and the heavy charged particles will be observed in the muon system. In 



the alternative DFSZ axion models [pq| , the stau decay length can be as short as 0(1) m 2£]. 
In such a case, the stau decay into an axino and a tau lepton may leave a kink signature. 
Hereafter, we do not discuss the axino LSP case for simplicity, keeping in mind that the 
models with the axino LSP in the DFSZ axion models can be studied in the same way as the 
gravitino LSP case. 

3.2 R-parity Violation and Cosmological Constraints 

If the R-parity is violated, the following terms are allowed in the superpotential: 

WRpV = -XijkLiLjEk + X'ij k LiQjDk + KiLiH u + -X'-^UiDjDk, (3) 

where we use the conventions, = —Xjik and X'-j k = —X'- k j. The first three terms violate 
the lepton number L, while the last term does the baryon number B. 
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Unless the lepton flavor violations are extremely suppressed, the B- or L-violation would 
erase the existing baryon asymmetry before the electroweak transition. Therefore, once we 
assume that the current baryon asymmetry is generated before the electroweak transition, 
the R-parity violations are severely constrained. The constraints are described as p], 
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E 



(0.6-l)xl0" 6 , (4) 
(3-6)xl0~ 7 , (5) 
(4-5)xl0~ 7 , (6) 



<(l-2)xl0- 6 (^) , (7) 



for the squark mass rriq ~ 200- 1200 GeV and the slepton mass mr~ 100-400 GeV. Here, fj, 
denotes the Higgsino mass parameter, and tan j3 is the ratio of the up- and down-type Higgs 
vacuum expectation values. The above severe constraints imply that, if R-parity is violated, 
the lightest superparticle among the superpartners of the Standard Model particles (which 
we call MSSM-LSP0) becomes long-lived and may leave interesting signatures at the LHC. 
We assume that the lighter stau is the MSSM-LSP and decays to Standard Model particles 
via the R-parity violating interactions. 

In our analysis the R-parity is considered to be violated only in the term ^\j^LiLjE^\ 
we assume the other sectors including the SUSY breaking terms have no R-parity violation. 
Then, the stau decays into 2- or 4-leptons. The decay channels for each type of the R-parity 
violating coupling are summarized as follows. 

(A sin 8)^ 

A123 : The decay rate of stau is r tot = m?, with Br(r — > ev) = Br(r — > uu) = 0.5. 

8-7T 

\ 2 n -\- gin 2 0) 1 
A^33 (i = 1,2) : The decay rate is r to t = 7773 m 7, with Br(r — > tu) = - — ■ — 2 



167T 1 + sin 

^ «in 2 ft 

and Br(r— >liu) 



1 + sin 2 ' 

\i3k (i = 1, 2; k = 1, 2) : Here is only one decay mode: r to t = r(r — > hv) = - — — — -m?. 
Ai2i;Ai22 : Stau decays into 4-body. 

Here, we ignore the mass of the leptons. A is the corresponding coupling constant, and v 



3 If gravitino is the LSP, the NLSP is the MSSM-LSP. If not, the LSP is the MSSM-LSP. 
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denotes both neutrinos and anti-neutrinos. 9 is the mixing angle of stau, which is defined as 

't\\ ( cos 6 sin#\ (t-C 



(8) 

t~2 J \— sm t> cos J \trJ ' 

where t\ and T2 are mass eigenstates, and t\ is now assumed to be the MSSM-LSP. Hereafter, 
for simplicity, we do not discuss the case where stau decays into 4-body, although it may be 
analyzed in a similar way as the 2-body case. 

For later convenience, we define an effective coupling A e ff, so that the decay rate becomes 

\2 

r tot = -s£-m=. (9) 
16vr r V ' 

That is, 

Ai23\/2sin6> for A123, 



A c ff 



A^fccos^ for A131, A132, A231, A232, (10) 

A i33 \A + sin 2 6 for Ai 33 , A 233 . 



Then the decay length cr of stau is estimated as 

CT = c (SfL m ~) ~ 990mm ( -^—) V m l \\ (11) 
V167T T ) V lxl ° / VIOOGeV/ v ' 

Therefore, if the coupling is A e ff < O(10~ 8 ), which satisfies the cosmological upper bound 
Eq. (U), the stau decay length is longer than 0(10) mm. 

Such small R-parity violation is cosmologically favored in the gravitino LSP sce- 
nario [30, 31 [ . When the gravitino mass is m^/ 2 > 5GeV, the R-parity violation of 



10~ 14 < Ajjfc, A' jfc < 10~ 7 naturally reconciles three well-motivated paradigms [31|: the 
primordial nucleosynthesis, gravitino dark matter, and the thermal leptogenesis. On the 
other hand, the relic gravitino decays via the R-parity violation, which may be observed by 



cosmic rays [30 — 32 1 . For m 3 / 2 > 100 GeV, the R-parity violation is constrained so tightly by 
cosmic ray measurements that the stau decay length becomes typically much longer than the 
detector size. For ?n 3 / 2 ~ 10 GeV, Ajj^ = 0(1O~ 8 ) can satisfy the cosmic-ray bounds, and 
the stau kinks can be observed. 

In our analysis we mainly concentrate on the LiLjE^ interaction for simplicity. When 
the R-parity is violated in the bilinear LiH n term, the stau decay causes a kink track, or 
splits into jet tracks, in the TRT, depending on the pattern of R-parity violation in the 
superpotential as well as in the soft terms. The former case can be described as the same 
picture as the Ajj^ case and thus covered by this paper. The signature of the latter case is a 
few jets plus 0-1 lepton, which is out of our scope and needs another method. 
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Table 1: The expected fractions of the daughter particles. The number ratios of e,/i and r 
produced at the r decay are given in the second group of the columns for each models. The 
particles/jet in the third group are observed in detectors, where the percentage shows the 
fraction of the number of each event. We use sin9 = 1 numerically for A133 and A233. 



Models 


e 


A* 


T 


e 


A* 


r-jet 


Gravitino LSP 








1 


18% 


17% 


65% 


A123 


1 


1 





50% 


50% 




Aj31 


1 








100% 






Aj32 





1 







100% 




Al33 


sin^ 9 





1 


59% 


9% 


32% 


A233 





sin 2 9 


1 


9% 


59% 


32% 



Table 2: mSUGRA parameters at the BC1 benchmark point. We will vary the R-parity 
violating couplings; since we consider very tiny R-parity violation, they do not affect the 
mass spectrum. 



M 


Afi/a 


tan j3 


Av 


sgn^t 


OGeV 


400 GeV 


13 


OGeV 


+ 



3.3 Model Discrimination 

Once the mother particle of the kink track is known, e.g. by studying decay chains of the 
SUSY events, the underlying model is shed light on by studying the daughter particle as 
well as measuring the lifetime and the mass spectrum. In particular, each model predicts 
a characteristic pattern of the branching ratio of the decay channels. Actually, when the 
R-parity is violated, the violation parameter determines decay products such as e, fj,, r and 
jets. When the stau decays into the gravitino, the tau leptons are produced. Since r decays 
into other leptons or hadrons, we expect to observe e, \x and jets at specific number ratios. 
Therefore, the reconstruction of the daughter particle is crucial for revealing the kink model. 
To be explicit, we summarize the branching ratios of the lepton production channels and the 
relative amounts of the measurable final-state particles for several models in Tab. [l]. 

4 Analysis 

We analyze the kink signature in the LHC by the Monte Carlo simulation. We consider the 
following three LHC setups: the center-of-mass energies are -Ecm = 7TeV, 8 TeV, and 14 TeV, 
with the integrated luminosity 2fb _1 , 5fb _1 and 10fb _1 , respectively. The high Pr one jet 
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Figure 2: The mass spectrum of our model point at the weak scale. 

plus large missing Ei trigger is used, which is standard for the SUSY searches. As described 
later, Pt > 120 GeV is required for the leading jet, and the missing is required to be 
larger than 100 GeV at the off-line analysis. These selections are related to the trigger and 
the trigger efficiency is higher than 90% for the events satisfying these conditions. The similar 
trigger condition will be used for the luminosity of 10 32 to 10 33 cm~ 2 s~ 1 for the high Pt one 
jet plus large missing E^ trigger. The single jet trigger is also useful for this analysis, and 
we can make use of the trigger redundancy. 

The benchmark point BC1 |33| is used with R-parity violating couplings varied The 
parameters are summarized in Tab. ^ and the mass spectrum at the weak scale is shown 
in Fig. §. The MSSM-LSP is the lighter stau with a mass of 140 GeV. Note that the tiny 
R-parity violation, one of our concern, hardly affects to the mass spectrum. The SUSY 
events are generated mainly via the squark-squark production channel and the squark-gluino 
channel; especially, the former dominates at the 7TeV LHC, while in the other setups they are 
in the same order. We utilize SUSY-HIT 1.3 |34|] to calculate the SUSY mass spectrum and 
the decay table. The SUSY events are generated by PYTHIA 6.4.23 |35| |, and the detector 
simulation relies on PGS4 Q. In addition, we use TAUDLA 2.9 |T?J (combined with TAUDLA 
C++ Interface |3^]) to simulate the decay of tau leptons, and PYTHIA to successive decays 
of the 7T° and r\ mesons. 

For the cases where the stau emits e or /x, the following selections are applied. 

4 Although the gravitino LSP scenario, one of our concern, is usually realized in the GMSB framework, we 
use this benchmark point in the analysis. The result is expected to be insensitive to the detail of the mass 
spectrum. 
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1. At least one jet with Pt > 120 GeV and the missing Et > 100 GeV are required. 
This off-line selection is related to the jet plus missing Ex trigger. The missing Et 
is calculated from the energies deposited on the calorimeters and the muon momenta. 
If the long-lived stau passes through the muon chamber and it is fast enough to be 
observed as a muon (/3 > 0.8), it is treated as a visible particle. 

2. \tj(t)\ < 0.63. 

3. P t (t) > 100 GeV. 

4. The stau should decay in the TRT, but before the 3rd module, i.e. \z\ < 712mm and 
563 mm < R < 863 mm. 

5. The kink angle 0.1 < A(j) < vr/2, where A(j> is the difference of the azimuthal angles 
between r and the charged daughter particle. 

6. The charged daughter particle (e or fi) should be reconstructed as a track. To this 
purpose, at least it should fly through the outermost layer of the barrel TRT (i.e. reaches 
R = 1066mm satisfying \z\ < 712mm). The reconstruction efficiency for it is set to be 
0.7 if P T > 20 GeV, or 0.6 if > 10 GeV; otherwise it fails this selection. 

When the daughter particle is a tau lepton, the above selection strategy should be modified 
slightly because the daughter tau decays promptly. The tau decays into leptons or hadrons, 
and we take both cases into account for the analysis. Especially, when the tau decays 
hadronically, it can generate multiple charged particles. It is noteworthy that, since the 
tau is highly boosted in the r decay, they fly into a relatively narrow conical region. Thus, 
we replace the last two items above by the following two conditions. 

5'. The "daughter cone" (see below) should be separated away from the track of the mother 
particle by 0.1 < Acft < ir/2 with respect to the azimuthal angle. 

6'. Among decay products of the tau, at least one charged particle should be reconstructed. 
The criterion is the same as the previous cut 6. 

Here a "daughter cone" is defined as the smallest cone which covers all the tracks of the 
charged daughter particles of the tau. Note that in the case of the 1-prong tau decay, the 
cone corresponds to the charged daughter particle, and the separation condition 5' is reduced 
to the kink angle condition 5. 

The numbers of the remaining SUSY events (or stau tracks in group II) are summarized 
in Tab. ||, in which decay length cr of stau is assumed to be 400 mm. This decay length 
corresponds to A e g = 1.3 x 10~ 8 for the R-parity violation models, and m 3 / 2 = 0.35 keV 
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Table 3: The cut flow for cr = 400 mm at the BC1 mass spectrum. The number of events is 
displayed in the group I, III and III', while the number of the stau is shown in II and II'. The 
daughter particle is supposed to be e//i in II and III. When the daughter particle is r, the last 
two items of II are replaced by II', and the result becomes III'. Note that one event contains 
two staus. The reconstruction efficiency of daughter track is taken into account. "Triggered" 
means the number after the cut 1, and the trigger efficiency (> 90%) is not included. 







7TeV,2fb- 1 


STeV^fb- 1 


14TeV, 10 fb- 1 


I 


total SUSY event 


673 events 


2832 events 


42463 events 


triggered event (cut 1) 


426 events 


1938 events 


36240 events 


II 


r track 


852 


3876 


72480 


r?(r) < 0.63 (cut 2) 


409 


1748 


28535 


Pr(f) > 100 GeV (cut 3) 


378 


1641 


26642 


f decay in TRT lst/2nd (cut 4) 


67 


230 


3642 


kink 0.1 < A<p < vr/2 (cut 5) 


46 


179 


2601 


daughter reconstructed (cut 6) 


28 


101 


1586 


III 


event with 1 or 2 kink 


24 events 


100 events 


1563 events 


event with 2 kinks 


4 events 


1 event 


23 events 


for the case where the stau emits r: 


II' 


separation 0.1 < A(p < ir/2 (cut 5') 


52 


189 


2805 


daughter reconstructed (cut 6') 


26 


95 


1391 


III' 


event with 1 or 2 kink 


24 events 


92 events 


1374 events 


event with 2 kinks 


2 events 


3 events 


17 events 



for the gravitino LSP models, for the stau mass = 140 GeV in the present model. (See 
Eqs. (§) and (|TT|).) The selections on the staus are employed in the group II (or II' when a 
tau lepton is emitted by the stau); the last numbers of the group II denote the numbers of 
the stau kinks at various center of mass energies. Finally in the group III (or III'), we show 
the number of the kink events, that is, events with at least one stau kink, and with two kinks. 

We find that, in this value of cr, several percent of the total SUSY events could be 
detected. In particular, the requirement that the stau should decay in the barrel TRT (the 
4th item of the group II) significantly reduces the events. We can enlarge the acceptance by 
including the staus decaying in the end-cap TRT, but not yet performed in this note. Even 
using only barrel TRT, the number of the kink events remains 0(1%) of the total SUSY 
events. 

As we discussed, when the daughter charged particle is a tau lepton, we slightly modify the 
selection strategy. However, to our delight, most of the daughter cones have small opening 
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Table 4: The superparticle masses at the weak scale and the total SUSY cross section for the 
model points considered in Fig. ||. We vary only M 1/ / 2 , and other parameters are the same as 
Tab. |. 



M 1/2 
(GeV) 


Masses (GeV) 


Cross section (fb) 


r 


9 


8TeV 


14TeV 


300 


103 


715 


2.95 x 10 3 


1.81 x 10 4 


400 


140 


932 


556 


4.22 x 10 a 


500 


176 


1145 


143 


1.31 x 10 3 


600 


212 


1355 


44.5 


472 


700 


248 


1562 


17.8 


194 


800 


283 


1768 


6.12 


87.1 



angles because of the boost; it was checked that the peak of the azimuthal opening angle 
distribution sits around 0.01 rad and almost all are < 0.1 rad. Therefore, the result appears 
very similar to the case where the daughter leptons are e or fi. 

In Fig. [|, we present the dependence of the number of the kink events on the stau decay 
length cr for 8TeV, Sib" 1 setup. Here, the daughter particle is assumed to be an electron, 
but the result does not change if it is a muon or a tau lepton. The figure tells that the number 
greatly depends on the decay length. This is because we demand that staus should decay in 
a very restricted region. 

In the figure, we also changed the SUSY parameter M]/ 2 . We show the masses of the 
relevant superparticles and the SUSY total cross section in Tab. |j. From Fig. |3] and Tab. ||, 
we can see that the ratio of the final number of the events to the number of the total events, 
for a given cr, is independent of Mi/ 2 - In fact, the number of the kink events is mainly 
controlled by the total SUSY cross section and the decay length cr. Here we can see that, 
when the stau mass is 100-200 GeV, the kink events can be discovered in a wide region of 
the decay length. 

We also show the same result for the 14TeV, 10 fb" 1 LHC setup in Fig. [|. The cr 
dependence is quite similar to the 7TeV case, but the number significantly increases due to 
the large cross section. If the decay length is O(10- 10 5 ) mm, we can achieve the discovery 
of the kinks in a wide region of the SUSY parameters. 

In Tab. ||, the accessible range of the lifetime are shown for the setup 8TeV 5fb _1 and 
14TeV 10 fb" 1 . We have also shown the corresponding ranges of the gravitino mass ?n 3 / 2 
for the gravitino LSP scenario, and that of the effective R-parity violating coupling A e g 
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Figure 3: Dependence of the number of the kink events on the stau decay length cr for 
8TeV, 5 ft) -1 . We also varied M 1/2 , one of the parameters of the mSUGRA model; the other 
parameters are the same as Tab. ||. The daughter particle is assumed to be an electron, 
though the result does not change if it is a muon or a tau lepton. In the simulation we 
generate 200000 events for each point to reduce the statistical uncertainty. 




Figure 4: The same as Fig. ||| but for 14TeV, 10 fb 1 . 
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for the models with R-parity violation. For 14TeV 10 fb 1 (8TeV 5fb 1 ), the ranges are 
0.03 keV < m 3/2 < 8keV (0.04 keV < m 3/2 < 2keV) and 4 x 10~ 10 < A cff < 8 x 10~ 8 
(1 x 10~ 9 < A cfr < 6 x 10~ 8 ), respectively. 



5 Discussion/Conclusion 

In this paper, we have investigated the possibility of probing SUSY models by analyzing the 
kink signatures of charged tracks in the TRT of the ATLAS detector. As underlying models, 
we discussed two models; (a) SUSY models with a gravitino LSP and a stau NLSP, and (b) 
R-parity violating SUSY models with a stau (N)LSP. It was shown that, if the decay length 
is 0(10- 10 5 ) mm, a large number of kink events can be discovered in a wide range of the 
SUSY parameters. For relatively light SUSY particles, i.e., gluino mass < ITeV and stau 
mass < 200 GeV, kinks can be observed already in the early stage of the LHC. 

In the case of gravitino LSP scenario, the kink signature can be observed for 
O(0.01 keV-lkeV) (see Tab. |). As discussed in Sec. [O], this corresponds to the most 



challenging gravitino mass range from cosmological point of view || . The reheating 
temperature after inflation must be O(100 GeV) or lower unless there is a large entropy 
production. To put it the other way around, if such a kink signature is indeed observed and 
if it becomes likely that the underlying theory is a gauge-mediated SUSY breaking model 
with the gravitino LSP, it has a significant impact on the cosmology. In particular, both of 
the inflation and the baryogenesis should occur at a very low energy scale. 

When the R-parity violation is a source of the kink track generation, the trilinear coupling 
is in the range of O(10" 10 -10~ 8 ). Although the MSSM-LSP can not be a dark matter, the 
gravitino LSP is a viable dark matter candidate. In this case, the severe constraint from 
the big-bang nucleosynthesis due to the superparticle decays can be avoided because of the 
R-parity violation. If the gravitino has a mass of ~ 10 GeV, the thermal leptogenesis can 
work successfully to generate the baryon asymmetry of the universe without suffering from 
the baryon washout. The decay of the gravitino might be detected by measurements of the 
cosmic rays, depending on the gravitino mass as well as the pattern of the R-parity violating 
parameters. 

Once the kink signature is discovered, the next step will be the discrimination of the kink 
models. Characteristics of the model are embedded in the final state particles of the kink 
decay. For instance, if we turn only on the coupling in the R-parity violating models, 
the final state is an electron (see Tab. |]). For ct = 400 mm, it is expected from Tab. ||] that 
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Table 5: The accessible range of decay length (more than 10 kinks can be observed) and the 
corresponding ranges of R-parity violating coupling and the gravitino mass for the models in 
Tab. |, at 8TeV 5 ft- 1 and 14TeV 10 fb -1 , 



M 1/2 
(GeV) 


Masses (GeV) 


N> 10 at 8TeV 5ft.- 1 


r 


9 


decay length (mm) 


m 3/2 (keV) 




300 


103 


715 


30 - 5 x 10 4 


0.04 - 2 


(0.1 - 6) x 10~ 8 


400 


140 


932 


50 - 1 x 10 4 


0.1 - 2 


(0.3 - 4) x 10" 8 


500 


176 


1145 


100 - 2 x 10 3 


0.3 2 


(0.5 - 2) x 10~ 8 




M 1/2 
(GeV) 


Masses (GeV) 


N > 10 at 14TeV 10 ft" 1 


T 


9 


decay length (mm) 


m 3/2 (keV) 


^cff 


300 


103 


715 


10 - 5 x 10 5 


0.03 - 6 


(0.04 - 8) x 10~ 8 


400 


140 


932 


20 - 2 x 10 5 


0.08 - 7 


(0.07 - 6) x 10~ 8 


500 


176 


1145 


30 - 6 x 10 4 


0.2 - 7 


(0.1 - 4) x 10~ 8 


600 


212 


1355 


40 - 2 x 10 4 


0.3 - 8 


(0.1 - 3) x 10~ 8 


700 


248 


1562 


60 - 1 x 10 4 


0.5 - 7 


(0.2 - 3) x 10~ 8 


800 


283 


1768 


90 - 4 x 10 a 


0.9 - 7 


(0.3 - 2) x 10~ 8 



we observe ~ 20 electron events for the LHC at 7TeV, 2ft- 1 , and ~ 100 for 8TeV, 5 ft- 1 at 
the BC1 benchmark mass spectrum. At the 14TeV run, the number increases up to ~ 1000 
for the integrated luminosity 10 ft- 1 . On the other hand, when we switch on Aj32 instead 
of Aj3i, the final state becomes a muon, which can be distinguished from the electron events 
by the detector signatures. Thus, it is expected that these models can be discriminated even 
in the early stage of the LHC run and much better for 14TeV. In the case of the tau final 
state, e.g. in the gravitino LSP case, ~ 20 r-jet events could be observed in the early LHC, 
while the numbers of the e and /i events are at most a few. In the 14TeV run, the number 
significantly increases, which enables us to distinguish the e, fi and r final-state models. 

In the case of the R-parity violating scenarios, our study covers the L^LjE^ interactions 
and generally, the bilinear L{H n terms. The LiQjD^ interactions, as well as some case with 
the bilinear violation, make staus decay into two jets plus or 2 lepton(s), which are not 
analyzed in this paper. The stau decaying hadronically in detectors is an interesting subject 
to study, which would complement our work. 
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